Energy gap of the bimodal two-dimensional Ising spin glass 
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An exact algorithm is used to compute the degeneracies of the excited states of the bimodal Ising 
spin glass in two dimensions. It is found that the specific heat at arbitrary low temperature is not a 
self- averaging quantity and has a distribution that is neither normal or lognormal. Nevertheless, it 
is possible to estimate the most likely value and this is found to scale as L^T~^ eyq>{—AJ/kT), for a 
Lx L lattice. Our analysis also explains, for the first time, why a correlation length ^ ~ exp(2 J/fcT) 
is consistent with an energy gap of 2 J. Our method allows us to obtain results for up to 10'' disorder 
realizations with L < 64. Distributions of second and third excitations are also shown. 
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In spite of its comparative simplicity, the two- 
dimensional bimodal short-range Ising spin glass model 
remains an interesting source of controversy. Although 
it is now widely accepted that the spin glass only exists 
at zero temperature [ij, l2|, the nature of excitations, in 
particular the energy gap with the first excited state, has 
commanded much interest in the literature. 

The bimodal model has bond (nearest-neighbor) inter- 
actions of fixed magnitude J and random sign. If we 
think of an infinite square lattice, without open bound- 
aries, it is easy to appreciate that any finite number of 
spin flips cannot result in an excitation energy of less 
than 4J. Nevertheless, some 20 years ago, Wang and 
Swendsen jSj] gave credible evidence that the energy gap 
with the first excited state should be 2 J in the thermo- 
dynamic limit. These excitations must involve an infinite 
number of spin flips. The issue here is the noncommuta- 
tivity of the zero-temperature and thermodynamic lim- 
its. In such a situation it is imperative to perform the 
thermodynamic limit first. 

Support for the 2 J energy gap has included studies 
that have involved exact computations of partition func- 
tions m, a worm algorithm [5| and Monte Carlo simu- 
lation [6[. Challenges to 2 J have also appeared. Saul 
and Kardar [7| maintained that the energy gap should 
be 4 J as naive analysis suggests. More recently [8|, it 
has been reported that the specific heat should follow a 
power law in temperature. In particular, it was proposed 
that the critical exponents must be the same as those for 
the model with a Gaussian distribution of bond interac- 
tions, indicating universality with respect to the type of 
disorder. Fisch [9|], using a steepest descent approxima- 
tion, has also suggested power law behaviour albeit with 
a different exponent. 

An important quantity, intimately involved with stud- 
ies of the thermodynamic limit at fixed temperature, is 
the correlation length. Essentially this measures the spa- 
tial extent of the influence of one spin on others. The 
correlation length is inflnite at a critical temperature. 
For the spin glass model of interest here, the correlation 
length has been determined by reliable Monte Carlo tech- 



niques [2,llO[ to be ^ ~ exp(2J/fcr). This is in agreement 
with Ref . [7[ and consistent with a qualitative study 11 1 . 
If this is true then hyperscaling predicts that the energy 
gap should be 4J. A 2J gap would be consistent with 
^ '^ exp{J/kT) as proposed in Ref. [g|. Another scenario 
is power law behaviour [8|. 

A good current review of the issues involved here has 
been given in Ref. [l2] ■ The power law behaviour of Ref. 
[8| is discussed in the light of new Monte Carlo data. 
The conclusion is that the suggested universality cannot 
be reliably proven with the computational facilities cur- 
rently available. Further, it is stated that extrapolation 
of the data of Ref. [f| to zero temperatures may not be 
plausible. The main message of this letter is to argue, for 
the flrst time, why a correlation length ^ ~ exp(2J//cT) 
is perfectly consistent with an energy gap 2 J, in apparent 
violation of hyperscaling. 

We have performed exact calculations of the degenera- 
cies of excited states at a fixed arbitrarily low tempera- 
ture. Each disorder realization consisted of a frustrated 
L X L patch with periodic boundary conditions in one 
dimension, embedded in an inflnite unfrustrated environ- 
ment in the second dimension. This choice of boundary 
conditions, with even L, deflnitely does not allow any 
first excitation with energy gap less than 4J. There are 
no open boundaries and no diluted bonds. A 2 J energy 
gap can only arise in the thermodynamic limit. 

Any planar Ising model is isomorphic to a system of 
interacting fermions. The Pfafflan formalism 11311 is par- 
ticularly convenient for disordered systems IJ, ll5J . Each 
bond is decorated with two fermions, one either side, so 
that a square plaquette has four; left, right, top and bot- 
tom. The partition function is given by Z ~ (detZ))^/^ 
where D is a real skew-symmetric 4A^ x 4N matrix for 
a lattice with N sites. Perturbation theory is used to 
determine ground state properties. Basically we require 
the low-temperature behaviour of the defect (meaning 
zero at zero temperature) eigenvalues of D. Exactly, 
D = Do + SDi where S = 1-t with t = tanh(J/fcr). 
The temperature dependence appears in 6 only. Dq is 
singular in the ground state and has eigenvectors | d) 



corresponding to zero eigenvalue, that is Dq \ d) = 0, lo- 
calized on each frustrated plaquette. These eigenvectors 
form the subbasis for the perturbation theory. 

At first order, the matrix Di, which is 2 x 2 block 
diagonal across bonds, is diagonalized in the subbasis. To 
continue at second order requires the continuum Green's 
fimction [ll,[3 Gc = (1 - P)gc(l - P) where P = Ed I 
d){d I and gc = gd +gc2- The first term gd is 4 x 4 block 
diagonal in the four fermions in each plaquette and allows 
us to connect frustrated plaquettes across two bonds. At 
second order we diagonalize D2 = DigdDi. The matrix 
5c2 is, just like Di, 2x2 block diagonal across bonds 
and only relevant for excited states. For higher orders 
we require Green's functions Gr constructed at previous 
orders. At third order the matrix to be diagonalized is 
Dj, = ^2(1 + GiDijgciDi and, generally at arbitrary 
order Dr = -Dr-i(l + 0^-2-0^-2) •••(! + GiDi)gciDi 
until all degeneracy is lifted. Here we define, for r > 1, 
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where Dr \ a^) — —e\ \ (31) and Dr \ PI) — e* | oi\.) with 
el ^ and there are N{r) such pairs at order r. Note 
that Gr is real here, as are all matrices; in contrast to 
Ref. [15| where they are imaginary. 
The internal energy of the system is 
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where cti is an Ising spin, C^ represents the sign of the 
bond and the nearest-neighbor correlation function can 
be written 14 1 
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where Q '" means the matrix element of ^ = —D~^ be- 
tween the two nodes decorating the bond {ij) \u\ . 

The full Green's function Q can be expanded for a finite 
system in powers of 5 
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where 



is the highest order of perturbation theory 



required. It is obvious that K,n, for to > 1, has no direct 
physical meaning. We provide here a brief outline of how 
the perturbation theory for the Green's function allows 
us to expand the internal energy exactly. 

Equating powers of 6 in QD = — 1, we obtain (for m < 

'^max ) 



DoK,n + DiK, 



m+l 



-OrnO- 



(5) 



Since PDq = 0, we can show that for to < 0, 

X™ = (1 + KiDi)gc{S,nO + DiK^+i). (6) 



It also can be proven that (l + KiDi) is idempotent, that 



{l + KiDif = il + KiDi). 



(7) 



Since both gc2 and Di are bond diagonal and (7c2-Di = 
i 0, we get Ko = {l + KiDi)gdil+DiKi)+g,2 + ^Ki, 
and for to. < 

K„,^{{l + KiDi)gciDi + ^)K^+i. (8) 

Substituting t = 1 — 5 and the Green's function in Eq. 
Q , the correlation function can be expressed in terms of 
Ki as 

00 
Cy(^.fT,) = (1 - 2K-+) ~ 2 Y, -^"[((1 + KiDi)gdDi 

m—l 

+ lr'-\l + K,D,)gd{l + D,K,)]-+ (9) 

We can now use the binomial theorem to expand in 
exp(— 2 J/fcT) instead of S, which eliminates the explicit 
effect of gc2, and the internal energy can be expressed as 



U=Y, e-"^"/'=^C/„ 



(10) 



m=0 



where Uq = -2JN + ^Jj^Uj) K^^ and for to > 



[/„ =2™+ij^[((l + Xii?i)g,ii?i)'"-^ 

[l + KiDi)gd{l + DiK^)]-+ . (11) 

Since D^^ — —1, it is obvious that 2E;--v A '" — 

-Efe)(A^)++ -E(.,)(^A)~ = -Tr(ADi) for any 
skew- symmetric matrix A. We can then show under the 
trace that for m > 0, 

(7™ = -2™ JTr[((l + Xii?i)g,iDi)'"(l + K^D^)]. (12) 

This can be rewritten using the idempotence relation 
(Eq. dlD) to get 

[/„ = -2™JTr(iPig,i(l + Z?iifi))™. (13) 

Finally, with recursive expansion, this can be expressed 



as 



[/,„ = -2"JTr(I?i.g,i(l + Dii^i))" 

= -2™ JTr(i?i.g,i(l + Z?iGi)(l + 152^^2))'^ 



ff^ TT"^ D^T- 



= -2"JTri?^ 
where 
R = Digdil + DiGi){l + D2G2) . . . (1 + £',,_G,„ 



(14) 



(15) 



This is sensible since it includes all the features of the 
ground-state calculation. It is also exact. 

We can imagine coloring plaquettes black and white; 
like a chess board. Matrices Dr and Gr are color diago- 
nal for even r: otherwise off-diagonal. Since gd is color 
diagonal, R is off-diagonal and it follows that t/,„ = 
for odd m. This explicitly excludes any 2J excitations 
involving a finite number of spins. 

The specific heat per spin can be derived in terms of 
the internal energy as 



1 dL/^ 1 ("14^ V^e-^^-A^C/. 
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We denote the degeneracy of the i excited state as Mi. 
Expanding \nZ, we obtain, for example, U2 = iJ{^), 
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As a means of establishing bearings, we first report 
a study of the first excitations of the fully frustrated 
Villain model [H, [l3|. Fig. [T] shows clearly that the 
-L^lnL. The correlation length is known 



ratio 



Ml 

Ma 



to be C ~ 5-^ 12Q|. Replacing InL with -\n5 - T^^, 
the correct form for the specific heat is obtained, that is 
Ci, ^ r~'^exp(— 4J/fcT). The degeneracy per spin of the 
first excited state is infinite in the thermodynamic limit, 
although only weakly (logarithmically) so. 
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FIG. 1: The degeneracy of the first excitations for the Villain 
model as a function of lattice size L. The curve is a guide for 
the eye. 

In Fig. [5] we show the distribution of ^ for the spin 
glass. It is clear that the most likely value scales as L^ . 
In this case we have an extra factor L, unlike InL for the 
Villain model. In consequence, taking ^ ~ exp(2 J/fcT), 
the specific heat varies like Cy ~ T~^ exp{—2J/kT) ex- 
plaining why a 2 J energy gap arises from the form of the 
correlation length obtained from Monte Carlo calcula- 
tions [2, llO|. Hyperscaling fails here. We also emphasize 



that the Villain model is not a spin glass and compar- 
isons, such as in Ref. [7|, are not meaningful. 

The distributions are scaled with a factor of L, prob- 
ably indicating that the amount of effort required for an 
experiment to find the mode of the distribution scales 
like L. This is consistent with the known difficulties that 
arise when trying to extrapolate Monte Carlo data to 
zero temperature. Also, the distributions are obviously 
neither normal or lognormal. Further, the data is not self 
averaging. In fact, the relative variance grows quickly 
with L which is unusually severe; convergence to a con- 
stant is normally expected 21 1 . 
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FIG. 2: (color online). Distributions, scaled with a factor L, 



for 



A/o U 



with various values of lattice size L. Each distribu- 



tion includes 10^ disorder realizations. 

We can discuss our choice of correlation length in the 
light of Ref. [ij]. Two crossover temperatures arc de- 
fined. Below T*{L) the ground state behaviour domi- 
nates. Above the finite-size crossover temperature T^{L) 
no size limitations are expected. It is clearly stated that, 
in all situations, T^ > T*. This immediately rules out 
^ r^ eKp(J/kT). For our case, we can in fact make the 
definite prediction T^ = 1.5T*. A power law behaviour 
for ^ must also be ruled out if we also expect similar 
behaviour for the specific heat. 

Pfaffians are also computed, for the partition function, 
in Ref. [J| at definite finite temperatures. For compar- 
ison, we have computed the mean of ln{jj^) and find 
fairly close agreement. Nevertheless, we do not believe 
that this is physically meaningful. The distributions of 
^ are not lognormal and the most likely value scales as 
L^, not L'*. We emphasize here that we are not comput- 
ing the entire Pfaffians. 

We have also checked distributions for second excita- 
tions. The most likely value of |^ — N scales as L^. 
We have subtracted off here the (infinite) total number 
of spins N. Single spin flips can occur anywhere and it 
is sensible to measure the internal energy relative to the 
unfrustrated system. However, the appropriate contribu- 
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FIG. 3: (color online). Distributions, scaled with a factor L, 
for the second term in the specific heat. Each includes 10^ 
disorder realizations, except for L = 64 which has 35, 000. 



state. All other treatments, except Ref. [7|, have re- 
quired extrapolation from definite finite temperatures. 
We have argued that an energy gap of 2 J is consistent 
with f ^ exp{2J/kT) as found from Monte Carlo sim- 
ulations. The manner in which our model is arranged 
excludes the possibility of any 2 J excitation for a finite 
system. Nevertheless, it should also be interesting to 
study systems with obvious 2J excitations, for example 
the hexagonal lattice or square lattice with open bound- 
aries or diluted bonds, and investigate the distributions. 
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FIG. 4: (color online). Distributions, scaled with a factor 
L, for the third term in the specific heat. Each includes 10^ 
disorder realizations. 
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tion to the internal energy is ^ 

in Fig. [HI the most likely value is close to zero. We note 
that the mode of lllf^-p')^ is roughly 0.002; indicating 
a close cancellation. Note also that the shapes of the 
distributions show clearly the dominant effect of the first 
excitations. The distributions for third excitations, that 
is of ^ - ^^ + KH)' - 2^: also show a most likely 
value close to zero as shown in Fig. [H Here, the mode 
of |j2. _ 2N scales as L^ and the unfrustrated system 
has M3 = 2A^. It is unlikely that higher excitations will 
interfer with our arguments based on first excitations. 

In conclusion, we have reported exact results for the ex- 
citations of the bimodal two-dimensional Ising spin glass 
by expanding in arbitrary temperature from the ground 
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